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Abstract 

We present a general method for studying long time asymptotics of nonlin- 
ear parabolic partial differential equations. The method does not rely on a priori 
estimates such as the maximum principle. It applies to systems of coupled equa- 
tions, to boundary conditions at infinity creating a front, and to higher (possibly 
fractional) differential linear terms. We present in detail the analysis for nonlin- 
ear diffusion-type equations with initial data falling off at infinity and also for 
data interpolating between two different stationary solutions at infinity. In an ac- 
companying paper |p, the method is applied to systems of equations where some 
variables are "slaved", such as the complex Ginzburg-Landau equation. 

1. Introduction 



The time evolution of many physical quantities is described by nonlinear, parabolic, 
partial differential equations. For most of these equations, to obtain a closed form 
solution seems to be a hopeless task. Therefore, one tries to determine certain qualitative 
properties of the solution, such as its existence and regularity for all times, or its long- 
time asymptotics. It turns out that, for certain equations, the long-time behaviour 
can be predicted because the solution becomes asymptotically scale-invariant; a trivial 
example is given by the usual heat equation: u = u" . The fundamental solution i) = 
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(47rt)~2e~~ is called scale invariant because its value at time t' may be obtained from 
its value at time t by scaling x and u. Moreover, for (suitable) initial data Mo(a;), 
the solution is asymptotically proportional to the fundamental solution, and therefore 
decays like . A less trivial remark is that the same asymptotic behaviour governs 
other equations; consider, for example, the heat equation with absorption: ii = u" — vP. 
Then, for p > 3, the long time asymptotics, given integrable initial data, is again given by 
the fundamental solution of the heat equation |T^ [T^. For long times, the nonlinear 
term vP has almost no effect. 

A similar situation is encountered in statistical mechanics in the theory of critical 
phenomena. Here too, it seems difficult to obtain an exact solution for many systems 
of interest. However, at the critical point, correlation functions become asymptotically 
scale- invariant; moreover, universality holds i.e. one may partition the set of statistical 
systems into different classes in such a way that, within a class, different systems have 
the same asymptotic behaviour at the critical point. 

In the theory of critical phenomena, scale invariance and universality can both be 
understood on the basis of renormalization group transformations. The idea is, very 
roughly, as follows: start from, say, a lattice system (like the Ising model), defined by 
its Hamiltonian H, and fix a parameter L. Then, perform the statistical sum over all 
degrees of freedom corresponding to fiuctuations of scale less than L. Rescale every- 
thing by L~^; we obtain a new lattice system defined by a renormalised Hamiltonian 
H' = Rl{H). This defines the renormalization group transformation. Now, iterate this 
procedure. The action of this transformation on the space of Hamiltonians is expected 
to have the following structure (this can be proven in certain cases): if we start with a 
Hamiltonian H which correspond to a critical point, the renormalization group will drive 
it towards a fixed point of that transformation. The fixed point is scale- invariant, which 
accounts for the asymptotic scale invariance of the original system, and a universality 
class corresponds to the basin of attraction of a given fixed point. Moreover, it often 
suffices to do a linear stability analysis of the RG transformation around the fixed point 
in order to determine qualitatively its basin of attraction. 

From a more technical point of view, the renormalization group allows the use of 
perturbation theory even where it seems to fail: in spin systems, there exist high and 
low temperature expansions which give good results away from the critical point but fail 
to converge up to that point. However such perturbative methods can be used to study 
the RG map itself, because the latter involves only the summation over a finite number 
of degrees of freedom (those of scale less than L). 

In 1978, Barenblatt wrote a book [Q] devoted to the role of scaling (which he calls 
"self-similarity of the second kind") in the asymptotic behaviour of solutions of PDE's. 
He suggested to look at these equations using the ideas being developed in field theory 
and in the theory of phase transitions. However, it is only recently, that, in several very 
interesting papers [l^, [l^, |l^ , it was shown how the RG idea can be used to study the 
long-time asymptotics of parabolic PDE's. The purpose of this paper is to extend this 
analysis and to show that it yields new mathematical results. For example, we can show 
that, for a large class of perturbations of the heat equation and for small initial data, 
the long time behaviour of the solution is given by the fundamental solution of the heat 
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equation. This result, which is related to those of [|TT|, Il2l IT6|, p[, also holds if we replace 



the second derivative in the heat equation by some other derivatives. 

We can also study what happens when one takes initial data that have a non trivial 
behaviour at infinity, i.e. such that uo{x) tends to two different stationary solutions of 
the PDE when x tends to infinity . We show again in a general situation, with small 
initial data, that the solution acquires in the long time limit a universal form, a universal 
"front", whose shape depends only on the boundary conditions and on the "universality 
class" of the equation. In an accompanying paper [^] we extend this result to the "slaving 
principle", i.e. to the fact that some (fast) variables in coupled systems are essentially 
following the behaviour of other (slow) variables. Those fast variables are an example 
of what is called relevant variables in RG terminology. As an example we consider in ||^ 
the complex Ginzburg-Landau equation (see also 0). 

The RG transformation for the PDE's is simply the integration of the equation up 
to an L-dependent time, followed by rescaling. After rescaling, we get a new equation. 
Upon iteration, this tends to stabilise to a fixed equation, whose scale invariant solu- 
tions are the fixed points of the RG transformation. Here also, the power of perturbation 
theory is enhanced by the RG approach: while the standard problem is to justify per- 
turbation theory for an infinite time, here we use it to study the RG map, which is only 
a finite time problem. However, as in the theory of critical phenomena, the RG idea is 
not limited in principle to the perturbative regime. 

We hope that the RG method will be a fruitful way to analyse a general class of 
equations. In any case, the equations discussed below provide very simple illustrations 
of rigorous applications of the RG method. We have tried to stress the analogy be- 
tween the partial differential equations analysed here and technically more complicated 
applications of the RG like field theory (see Section 5). 



2. The Renormalization Group Map 

We will explain now the RG method for an equation of the type 

u = u" + F{u,u',u") (1) 

(for notational simplicity we take the spatial variable x G R: everything below can as 
well be done in R"). We are interested in the asymptotics of the solution of the form 

M(2;,t) ~ rt/(xr^) (2) 

as t — s> cxD. 

The standard way to study this problem is, 12, 16], to first look for a scale- 
invariant solution (if F is scale-invariant) which reduces to solving an ordinary differential 
equation (see Section 5) and then to establish the stability of this solution. Usually, one 
uses a priori estimates like the maximum principle to prove stability. 

The RG method transforms the problem of large time limit into an iteration of a 
fixed time problem followed by a scaling transformation. The scale-invariant solution 
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emerges then as a fixed point of a map in the space of initial data, the RG map, and 
the stabihty analysis becomes the analysis of the stability of the fixed point under the 
RG. As wc will see, the method does not depend on any a priori positivity properties 
and is applicable to systems of equations and to equations whose leading term has other 
derivatives. 

We now explain the idea without spelling out any concrete assumptions on the F 
in (1) nor on the spaces and norms. The choice of the latter, as usual in non linear 
problems, depends on the particular problem. Thus, let us fix some (Banach) space of 
initial data S. It will be convenient to take the initial time to be t = 1. Next, we pick 
a number L > 1 and set 

UL{x,t) ^ L'^u{Lx,LH) (3) 

where a will be chosen later and u solves (1) with the initial data f E S. The RG map 
R: S ^ S (this has to be proven!) is then 

{Rf){x) = UL{x,l). (4) 

Note that ul satisfies 

iiL = ul + Fl{ul,u'l,uI). (5) 

with FL{a, b, c) = L'^+''F{L-''a, L-^-^b, L'^-'^c). 

We may now iterate R to study the asymptotics of (1). R depends, besides on a, on 
L and F. Let us denote this by Rl,f- We have then the "semigroup property" 

Rl^^f = Rl,f^„-i o ■ ■ • o Rl,Fl ° Rl,f- (6) 

Each R on the RHS involves a solution of a fixed time problem and the long time problem 
on the LHS is reduced to an iteration of these. Letting t — L^", we have 

u{x,t) ^t-^{RLn,Ff){xt--^). (7) 

Now one tries to show that there exists an a such that 

Fin F* , RL^^pf f* (8) 

where 

i?L,F*r - r (9) 

is the fixed point of the RG, corresponding to the scale-invariant equation u — u" + F*. 
Then, rescaling x, the asymptotics of the original problem is given by 

u{xt^^^,t) ^t-^f*{x). (10) 

We will now illustrate the RG method in two concrete cases. 
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3. Gaussian Fixed Point: Diffusive Repair 



For F = the equation (2.1) is just the diffusion equation and the corresponding 
fixed points are trivial to write down. We concentrate on integrable initial data; for a 
discussion of fixed points relevant to other data see the end of this section. 

Going to Fourier transform and putting a = 1, we have 

Rf{k) = e-''^'-'^-'U{L~'k) ^ Ihfik). (1) 
-Ro has a line of fixed points, namely the multiples of /q, with 

foik) = e-'\ (2) 

/g is of course the initial data of the scale invariant solution u{x, t) = (47rt)~2e~ 4t . 

Note that Ro has nice contractive properties, provided / has some smoothness. In- 
deed, let f = Jq + g with g{0) = 0. Then, Rof = /o + Rog and, if g is C^, we have 
|9'(x)l — ■'i^I^'(t)I' -^^^ < s < 1, which leads to contraction in many "natural" norms. 
Here, we shall consider the Banach space B of functions / with / e C^(R), equiped 
with the norm (we are inspired here by [Q) 

11/11 = sup(l + fc^)(|/>)| + |/'(fc)|). (3) 

k 

The 1 + factor provides smoothness for / (actually we only need 1 + with p > 3) 
whereas the /' term gives decay at infinity (e.g. xf{x) e L^). Now, for ^'(0) = as 
above, we get 

\\Rog\\<CL-'\\g\\, (4) 

with C independent of L. Here and below, C will denote a generic constant, which 
may change from place to place, even in the same equation. However, we shall write 
Cl;Cf,... for (generic) constants that depend on the choice of L, F etc. In the proofs, 
L will be chosen large enough so that e.g. CL~^ < 1. Then other quantities like e in 
Theorem 1 below are chosen small enough in relation with L. 

Next, we discuss the domain of attraction of the Gaussian fixed point (2) in the 
set of equations of the form (2.1). We will take the nonlinear term F to be a function 
F : — > C which is analytic in a neighbourhood of the origin. Note that for a monomial 
F(a, b, c) = a^'^b^'^c''^ we have Fl = L''^'' F with 

dp = ni + 2n2 + Sn^ — 3. (5) 

In general, we define dp for an analytic F by taking the smallest of the numbers (5) 
computed for the monomials in the Taylor series of F at with non- vanishing coefficients. 
We call F irrelevant, if dp > 0, marginal, if dp = and relevant, if dp < 0. 

A. The irrelevant case 

We can now state the main result for an irrelevant F. 
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Theorem 1. Let F : — > C &e analytic in a neighbourhood of with dp > and fix 
a S > 0. Then there exists an e > such that if \\f\\ < e, the equation 

u = u" + F{u, u', u") u{x, 1) = f{x) (6) 

has a unique solution which satisfies, for some number A — A{f, F), 

\mit'~'\\u{-tKt)-At-'^f*{-)\\=0. (7) 

t— >oo 

Remark 1. As will be evident from the proof, we could as well consider a more general 
class of equations 

ii^-i-Af/^u + F , k:R^^C (8) 

where {—A)'^/'^f{k) = \k\^f{k) and F is analytic in u and its partial derivatives up to 
order /3. f*{k) is then e"''^'" and u{x,t) ~ At~'^ f*{xt~'^), provided dp > where, for a 
monomial F — IljLi n'=i(5"^w)"''^ , we set 

j=i i=i 

We restrict < p. In (2.3) we set a — N and replace L'^t by L^t. In (3), we replace 
A;^ by \k\P with p> (3 + N. In (5) above, we had iV = 1, /3 = 2 and = 0, 1, 2. Of course 
we could also take in F fractional derivatives, \k\'^ in momentum space for 7 < 

Remcirk 2. The statement (7) translates in momentum space into 

\uik,t)-Ae-"''\ < Ct-'^+\l + t''kY\ (9) 

Taking smoother initial data improves the fall off on the RHS accordingly. In x space, 
the convergence in the norm (3) implies convergence both in and in L^. For L^, we 
use Schwartz' inequality to get 

\\f\\i<C\\{l+\x\)fh 

and Plancherel to bound 

II/II2 + Ik/lb <C||/|| 
For L-,weuse ||/||oo < ||/||i < C||/||. 

Proof. We start by discussing the local existence of the solution. Thus, turn (6) into 
an integral equation 

u{t) = e(*-^)^V + dse'^'F{t - s) = Uf + N{u) (10) 

J 
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where F{t) = F(u(t),u'{t),u"{t)) and d = We solve (10) using the contraction 
mapping principle by introducing the Banach space of functions u{x, t), t ^ [1, L^] with 
the norm 

= sup \\u(t)\\. (11) 

tG[l,L2] 

We shall show that T{u) = Uf + N{u) maps the ball Bf = {u\\\u — < 11/11} iiito 
itself and is a contraction there, for ||/|| < e = e{F, L) small. 

We need to estimate ||A^(u)||l. In this the use of momentum space in the norm (3) 
is very convenient: we are indebted to for this observation. Thus, Taylor expand F 
in (10) and take the Fourier transform: 

F = anU ^ *u' * u" (12) 

neN3 

where * is the convolution. Inserting (12) into (10) and using the bound on the deriva- 
tives \v^\k)\ < C\\u\\l{1 + k"^)-^ , for / = 0,1,2, we get 

dse-'^'\{u*''' * u'*"' * ^?™')(A:)| < (C||m|U)''"(- " ^ ^* ^ 



2 



k 

■Jil + ik- piY)-' . . . (1 + {pmyr'dpi ...dp^< 

CL[C\\uhni + k'')-\ (13) 

with m = rii + 71,2 + 71-3. The same estimate holds for the derivative with respect to k and 
we get a convergent bound on ||A^(ii)||L: since F is analytic, |an| < (Cj?)™, and hence 

\\N{u)h<CL,FMl (14) 

{dp > implies m > 2), provided is small enough. This holds in Bf, if ||/|| < 

e{F,L), since < C\\f\\ and so < (C* + 1)||/||- Thus T maps Bf into itself. 

In the same way, we estimate 

\\N{u,) - N{u2)\\l < CfL''{\\ui\\l + \\u2\\l)\\ui - U2\\l- (15) 
Thus, for ll/ll < e(F, L), (10) has a unique solution in Bf 

u{;L') = e^'^'-'^''f + v (16) 

with 



Now we write 



^11 < CLAfV (17) 



f = Aof*+go 
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with Aq — /(O). Note that /o (0) = 1, hence go{0) — 0. The reason for this decomposition 
is that Ro is a contraction when it acts on functions g with ^(0) = (see (4)). Also, 

ll^o|| = ||/-/(0)/oll<q|/||. (19) 

Then, 

Rf = Lu{L; L2) = Rof + Lv{L-) = AJ* + g, (20) 
where Rq is defined in (1), 

Ai = Ao + v{0) (21) 

and 

gi^R^go + LviL-)-viO)f;, (22) 
so that gi{0) = 0. We have then the estimates, using (17), 

l^-^l <Ci,F||/ir (23) 
\\Lv{L.)-v{0)fo\\<CL,F\\fr (24) 

and, combining with (4), we get 

\\g,\\ < CL-%o\\+C,M\fr < L-^'-'^\\f\\ (25) 

for ll/ll < e{L,F) and L large (depending on 6). 

The proof is now completed by iterating this procedure. Set 

f^ = RL„,Ff = Anf^+gn (26) 

So, with V above equal Vq, 

An+l = An + Vn{0) (27) 

gn+1 = Rogn + LVn{L-) - V^{Q) f* (28) 

and assume inductively that 

|A„|<(C-L-")||/|| , \\g4<CL-^'-'^^\\f\\, (29) 

so that ||/„|| < C||/||. Then, repeating the above analysis and noting that in (14) and 
(15), Cf is replaced by L-'^'^^Cf (see (2.5)), we get, instead of (17), (23), (25), 

\\vn\\<CL,FL-'''-\\f\\' (30) 

\A^+i-A^\<CL,FL-^'-\\ff (31) 
||^„+i|| < CL-i^„|| +Ci,^L-"'^-||/||^ (32) 
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Thus, (29) iterates {dp > 1) and we get that An A with \An-A\ < CL,FL-'"^^\\f\f < 
CL-""^^ 11/11 for 11/11 small, 

\A-f{0)\ = \A-Ao\<CL,F\\ff (33) 
and (see (2.7), (26)), for t = L^", 

\H-t'/^t) - Ar'^f*{.)\\ < ct-'^'^'Wfl (34) 

It is trivial to extend this bound to t = tL"^"-, with 1 < r < by replacing everywhere 
by rL^. This proves the claim. □ 

B. The marginal case 

Let us now consider the marginal cases, i.e. 

F = -u^ + G{u,u',u") (35) 



or 

F = 2uu' + H{u, u', u") = (m^)' + H{u, u', u") (36) 

where G and H are irrelevant. (36) is just the Burgers equation with a perturbation. 
Since the discussion of that equation uses some ideas introduced in the next section, we 
shall discuss it at the end of Section 4. In (35), we need the negative sign: for G = 
and a positive sign, any non zero initial data leads to a solution that blows up in a 



finite time [|18| , ^ . However, we do not have to assume that the initial value is pointwise 
positive; rather we want / to be near a small multiple of /*. Although a general analytic 
irrelevant G could be treated, we assume for simplicity that the Taylor expansion of G 
starts at degree 4 or higher. Without this assumption, in the proof below, we would 
need to consider first a "crossover" time during which the G term would dominate u^. 
We will therefore scale u by A2, where A is chosen so that u{-, 1)(0) = 1, and consider 
the equation 

ii = u" — Xv? + Gx{u, u , u") (37) 
with G\{z) = \^^G{\^ z) (which is of order X^/"^ for A small) and initial data 

/ = /o+^?, ^(0) = 0. (38) 

We have then 

Theorem 2. For any 5 > there exist Aq , e > such that for < A < Aq and \\g\\ < e, 
we have 

hm tH\ogty-'\\u{-t-2,t) - (-^tlogt)-5/*(.)|| = (39) 



t— >oo 



'2V3- 



TT 
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Proof. We proceed as in the proof of Theorem 1. We use (10) and we get, for A small, 
instead of (14), 

\\N{u)\\l<Cl,gX (40) 

{\\u\\l is of order one here). This, with the analogue of (15), proves local existence, 
leading to 

Rf ^ R^f + LviL.) 

with 

\M < Cl,g\ (41) 

Since is marginal, we have to handle rather explicitely its main effect on u which 
will lead to the logarithmic correction in (39). Let, in (10), N{u) = —XN^i^u) + Nc{u) 
corresponding to the two terms in (37). So, 

N3{u){t) = dse'^\u{t - s)f 
Jo 

We have 

\\Ng{u)\\l < Cl,gX^/^ (42) 

Denote by u*^ the solution of (37) with Gx = and / = Afg, and by ua the solution of 
(37) with/ = A/* + So, 

<(t) = Ae(*-^)^Vo-AA^3«)(i) (43) 

UA{t) = Ae(*-^)^Vo + e^'-'^'^g - XN^Mit) + NGMit) (44) 
Here A = 1, but we shall use below the following bound, for \A\ < 1 and \\g\\ < e: 

||iV3(^A) - iV3«)|U < CL{A^g\\ + \\gf) + Cl,gX^. (45) 
To prove (45), first show, like in (15), but with a cubic nonlinearity, that 

\\Ns{ua) - N3{u*a)\\l < ClUWaWl + \\u*A\\Lf\\uA - u*a\\l + Cl,gX^ 
Then, use (42), (43), (44) to show 

\\uA - u*a\\l < CMl + CL,GAi + A||Ar3(^x^) - Nsiu*A)\\L 

Finally, use \\uW\l < Cl\A\ and \\ua\\l < Cl{\A\ + \\g\\), which follows from (43), (44); 
(45) follows for A small. 
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If we define 

Jo 

we get, inserting (43) in N^{u*j^), and using the bound ||-/V3(w^)||l < Ci,|yl|^, 

\\N^{u*^){L^) - Ah*\\ < ClXA^ (46) 

for I A I < 1. Now, write 

V = -Xv* + w (47) 
where, using (42), (45), (46), for A = 1 and ||5f|| < e, 

\\w\\<CL,GX{\\g\\ + X^). (48) 

((46) gives a term of order here). Again, we write 

h = Rf^ AJ* + g, (49) 

with = 1 + v{Q) {Aq = 1 here, by (38)) and gi given by (22). Denote v*{Q) = /?. We 
get from (48) 

\A,-l + X(5\ = \wm < CL,GXi\\g\\ + X-^) (50) 

and, from (4), (41), 

\M < CL-'\\g\\ + CL,GX (51) 

The iteration is as follows: the G term will run down with a factor L^"'^'^, which 
improves the bound (42) and the last term in (45), but we shall see that, unlike the 
situation of Theorem 1, An will go to zero and we have to keep track of the correction 
of order A^ in (46). We use (26), (27), (28) and we have 

Vn = -XAlv* + Wn (52) 

with, using again (42), (45), (46), 

IknII < CL,GX{Al\\gJ + \\gnf + X^ L-^'^ + XA'J. (53) 

So, 

\An+i-An + XpAl\ < CL,GX{Al\\g4 + \\gnf + X^L-''''^ + XA^) (54) 

Using the fact that P > (see (58) below), we see that An decreases and so, using (28), 
(4), (52), (53), lit;* II < C, and Cl,gXAI < L-\ for A smaU, we get 

lbn+i|| < CL-%n\\+CL,GXiAl+\\gn\\' + X^L-^'<' + XAl). (55) 
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The iteration of An,gn, leads to 

42 - 



(56) 



ll^nll < CL,Gn-^ (57) 

for n > 1, with \bn+i — bn\ < CL,Gn^K which gives < Cl^gV^- Indeed, the leading 

term in (55) is and, in (54), Al\\gn\\ We get An = {2\pn)~^ + 0{n-^). Finally, 
we compute 

P=-^ ds [ dpdge(^-^^)(^^+''^+(^-^)^) = (58) 
47r2 7o J-R? 2^/2>^x 

Using 2nlogL = logt, the claim follows as in Theorem 1. □ 



4. Non- Gaussian fixed point: stability of a front 



If equation (2.1) has other (constant) stationary solutions than m = 0, we may study 
the problem where lim^^-too w(x, 1) takes two different values. This corresponds to a 
front in the initial data. One may then inquire about the stability of this front and 
the universal features of the long time solution in such a situation. These problems 
are discussed in general in [g, [7| and, for coupled equations related to the equations 
discussed here, see [|], |], ||. Here, the restriction to one space dimension is essential. 

An example of such an equation is the nonlinear diffusion equation 

u = diil + a{u,du))du) , uix,l) = (pix) , lim (j)ix) = u± (1) 

a;— > it 00 

where we denote d = This is a special case of (2.1), with F = au" + ^{u'Y + ^u' u" . 
If a is analytic and if a(0,0) = 0, F is irrelevant. (1) has u = const, as a stationary 
solution. 

Remark 1. The most general such equation in one variable we could deal with is 

u= {l + a{u,u',u"))u" + b{u,u',u"){u'f (2) 

with a and b analytic as before. See Remark 2 below. 

To understand the problem in the RG setup, let us first consider the trivial case 
a = 0. This is of course exactly soluble. We have 

u{Vtx,t + l) = / e'-^^y-'^^' (t){Vty)dy 

^t-*oo + {u+ - u^)e{x) = (pKx) (3) 
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12,-/ 

where e{x) — /^^ e'i^ In RG terminology, we have the "Gaussian" fixed point (f)^ 
corresponding to the u± boundary condition problem. It is easy to check that (pQ is a 
fixed point for the map 

RL(f>^uiL;L') , = (4) 

Note the absence of the multiplicative L factor in i?^ (in (2.3), we have a = 0): the 
initial data are not normalizable by their integral and u{x, t) does not, in general, go to 
zero as t goes to infinity. 

The stability of this fixed point is also easy to understand. We write = 0q + /, 
where /(±oo) = 0. Then v = v", with v{x, 1) = f{x), and the analysis of the previous 
section applies. Thus we expect 

u{x, t) ~ 0*(^) + ^/o (^) + 0{t'-') (5) 

This is the asymptotics we wish to prove for (1), except that (f)^ and /q are replaced by 
non trivial fixed points. We will prove a theorem for small data, so u± will be taken 
small. What wc do below is an (elementary) example of what is called, in the theory of 
critical phenomena, the "epsilon expansion". Before stating the precise results, let us 
do the heuristics. 

Upon scaling UL{x,t) — u{Lx, LH), ul satisfies (1) with a replaced by aL{u,u') — 
a{u, L~^u'). Thus, we will search for a fixed point for the RG (4) with u satisfying the 
equation obtained in the limit L — > oo: 

u^d({l + a*(u))du) , u{x,l)^(l)*(x) (6) 

with (p*{±oo) = u± and a*{u) = a{u,0). We will find the fixed point (p* by finding a 
scale invariant solution to (6): 



u{x,t)^cl>*{^). (7) 



We get (replacing by x) 



d{{l + a* {(/)*)) d(i)*) + ^xdcj)* = (8) 
and we look for a solution 

(/.* = <^*+^ (9) 
with ■0(±oo) = and 0q is the Gaussian solution (3). Thus, we get for the equation 

AV' = 9(a*((/)* + V)9(0* + V')) (10) 

with 

A^-d^ - ^xd. (11) 
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Let us also denote 

i;o = d{a*{cl>l)dcl>l). (12) 

We shall discuss the properties of later. We solve the fixed point equation (10) in 
the space of functions equipped with the norm 

U\\n= max sup|a'"V(a;)|e^. (13) 

0<m<N X 

Let the degree of a* (i.e. a*{u) = 0{u'^) for u small) be c? > 0. Then we have the 

Proposition. Let a : ^ C in (1) be analytic in a neighbourhood of and let N be 
an integer. Then there exists an e > such that, for \u±\ < e in {3) , (10) has a unique 
solution with 

\\^-ipo\\N<Ce^''^'. 



The proof of the proposition will be given after the one of Theorem 3 below. We will 
also see that ■0o can be written down exphcitely and is (9(e''+^) in the norm (13). 

Given (f)*, we return to equation (1), and write 

u{x, t) = 0*(^) + v{x, t) = u*{x, t) + v{x, t). (14) 

Because of (6, 7), v satisfies the equation 

v^d{{l + a)dv + {a-a*{u*))dvr) , v(x,l) ^ f(x) (15) 

with lima-^-i-oo f{x) = 0. Now the RG is as in Section 3: 

RLf^Lv{L;L^) = VL{;l) (16) 

However, the fixed point will not be the Gaussian /q of that Section, but a different one. 
We find it by taking the scaling limit: vl satisfies the equation 

vl = a((l + aL)dvL + L{aL - a*{u*))du*) (17) 

with aL = a{u* + L~^vl, L~^d{u* + L~^vi,)), and therefore we look for a fixed point 
satisfying the linear equation obtained from (17) by letting L — > oo: 

V* = d{d{{l + a*{u*))v*) + —{u*, 0){du*f). (18) 



This is solved by setting 



v*{x,t)^^^n^) (19) 
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whence, ^(^((l + a*{(j)*))f*) + \xf* + /i) = where h = 0(0*, O)(a0*)2. We solve 
this equation by integrating once (we shall look for an integrable /*, so the constant of 
integration is zero) and then solving a first order equation: 

fX 1 rx y + 2aa* {<t>* (y)) j 

f*(x) = { h{y)F{y)dy + J\f)e~^Jo . (20) 



I j-x y + 2da*{<j>*{y)) , 

where F{x) = —(1 + a*{(j)*{x))~^e ^ ■'o '^r^^WIW ^ , For future convenience, we normal- 
ize this by choosing M such that J f* = /*(0) = 1. Note that /* is, for 0* small, a small 
perturbation of /q. In particular, /* is smooth and decays rapidly at infinity. 

Remark 2. This is the only place where (1) is simpler than (2). For (2), (18) is replaced 

by 

(9(7* 

V* = {1 + a*{u*))d^v* + 2b*{u*)du*dv* + {—{u*){d\*) + 

ou* 

— («*)(a«*)V + ( — (^x*,0,0)(9V) + — («*,0,0)(a«*)2)(9^z*) (21) 

with a*{u) = 0(^,0,0) and b* similarily. v* will not be as explicit as (19,20) and ^ 
in front of /* in (19) will be replaced by since, while / v* is constant for the 



solution of (18) (integrate both sides), it is not conserved by (21); see |jT3|, |Tj] for a 
discussion of a similar effect. 

Before stating the main result of this section, we need to specify the space of initial 
data (j) of (1), i.e. the / of (15). We can not directly use the norms (3.3) introduced in 
Section 3. The reason is that the function (p* in (9) involves the Gaussian fixed point 
(3) which has no falloff at infinity. In particular it is not in the Banach space we used 
before. However, the only way u* enters in (15) is as du* that falls off like a gaussian, 
or by multiplying functions that we expect to fall off at infinity. As a consequence of 
the latter possibility, we may not use pure momentum space bounds as in (3.13): the 
convolutions would then involve u* which is too singular (|0q| ~ |^|~^ for small k). 
Instead, we introduce a norm that encodes more sharply than (3.3) both the long and 
short distance properties of the solution. 

Thus, let X be a non negative function on R with compact support on the interval 
(— 1, 1), such that its translates by Z, Xn = x(- — n), form a partition of unity on R. For 
/ G C^, we then introduce the norm 

11/11= sup (l + n^)(l + fc2)|x:3^/(fc)|. (22) 

neZ,feGR,i<2 

Roughly, 11/11 < 00 means that / falls off at least as at infinity and f{k) as /c~^. 
Note that the derivatives of 0* and /* and its derivatives have a finite norm. Comparing 
with (3.3), we have k'^ instead of k"^ but two derivatives act on /. We used k"^ to do 
the convolutions in (3.13). The could be changed to anything increasing not faster 

than e « (coming from (13)). These norms for different choices of x may be shown to 
be equivalent. 
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Theorem 3. Let a : — > C &e analytic in a neighbourhood of with positive degree 
and fix a 5 > 0. There is an e > such that for \u±\, \\f\\ < e, the equation 

ii = d{l + a{u,du)du) , u{x,l) = (p*{x) + f{x) (23) 

has a unique solution, satisfying 

lim t'-'\\u{Vt;t) - ct>*i-) - ^n-)\\ = (24) 
where f* is given in (20) and (jf in the Proposition. 

Remark 3. The convergence in the norm (22) again imphes convergence in and 
in see equations (39) and (40) below, apphed to i = 0. Equation (24) means 
that u{x, t) behaves, for |x| < C and t oo as (/)*(0) + At^^ + 0{t^-^), where A = 
(0*)'(O) + /(0)/*(0). Thus, locally, the solution goes to a constant i.e. a stationary- 
solution of (1), with a diffusive correction. This solution is "selected" by the boundary 
condition at infinity, u±, and by a. 

Remark 4. Since in (9) satisfies < Ce^ (this follows easily from the Proposition) 
we can replace 0* in (23) by 0q, given by (3), which makes the hypothesis more explicit. 

Proof. We consider the equation (15) and the RG (16). It is convenient to separate the 
/* piece from /. We write 



and correspondingly 



with 



f-mf* + g (25) 



v* + w (26) 



v*{x,t) = l^f*{^^). (27) 

First, we have |/(0)| < EneZ 1x^(0) | < EneZ iK? ^ ^WfW H/H ^ ^' ^^^^^ 
implies 

ll^ll < C\\f\\ < Ce. (28) 
The equation we finally study is the one satisfied by w: 

w = d'^w + K , w{x,l) = g{x) (29) 
where ^(0) = (see (25) and the normalisation in (20)) and where (see (15), (18)) 

K = d[(a - a*{u*))dv* + adw + {a- a*{u*) - -^{u*, 0){du*)))du* - v*da*{u*)] (30) 
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with a — a{u, u'), u — u* + v* + w. The exphcit form of K is not important for us, all 
we need to know is that it is given as a convergent power series in u*, v*, w and their 
derivatives up to second order, has no linear term and no term which is only a power of 
u*. We shall also use some scaling properties of K discussed below. 
In RG language, we need to control (see(16)) 

L"«;(L"x, L^") = RLn,K9 = i?L,K„_i ■ ■ ■ Rl,k9 (31) 

where now Kn+i is obtained from Kn by replacing (sec (16), (4)) 

d'v* L-'-'d'v* , d'w L-'-'d'w , d'u* L-'d'u*. (32) 

To see the behaviour of the terms in (30) under this scaling, define d — nQ + n^* + n^,, — 3 
where nj}*,nw count the number of v*,w factors and their derivatives, while ng is the 
total number of derivatives. We say that a term in (30) is irrelevant if <i > and marginal 
if d = 0. A term with d < would be relevant, but these terms have been included 
in the equations for u* and v*. To check this, notice that each term in (30) has either 
more than two derivatives or two derivatives and at least one v* or w factor, so that 
d>0. In fact, the only marginal terms in (30) are of the form {u*yd'^w or {u*ydu*dw. 
To check that there are no terms like that with w replaced by v, observe that such 
terms are cancelled because of the a — a* in front of dv* and because of the subtraction 
between the two last terms in (30). A term like d-^{u* ,0){du*y is also marginal, but 
it is cancelled by the subtraction in the factor multiplying du* in (30). 

Eventually, using ^(0) = 0, we will show that H-Rin^xfi'll goes to zero like L~^^''^^'^. 
Using (14), (26), (27), (31), this will prove (24). We split the proof again into three 
steps: the proof of local existence, the control of R and the iteration of R. 

(a) Local existence. We have to solve the integral equation 

w(t) = e^'-^^^'g + dse'^'Kit - s) 
Jo 

(with obvious notation) using the contraction mapping principle with the norm 

= sup ||w(-,t)||. (33) 

tG[l,L2] 

As in Theorem 1, we expand K as a power series and estimate a generic term 

a{x) = f~'ds I dye-''\x-y)u*{yyF{y) 

where F{y) is a (non empty) product of the functions v*, w and their derivatives, and 
derivatives of u*, and / > 0. We supressed the t — s variable in u* and F. We localize 
the y variable: 

«(^) = E ds I dye'^\x - y)Xra{y)u*{yy F{y) ^ ^ a^(x). (34) 

meZ meZ 
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We want to bound 

/3mn = sup + k^)Xnd'am{k)\. (35) 

k 

We distinguish between |n — m| > 2 and \n ~ m\ < 2. 

(A) Let first \m — n\ > 2. Then Xm and Xn have disjoint supports and e^^^{x — y) is 
smooth uniformly in s. We write 

Pmn = sup 1 1 dxe-'^'"'-' ds J dy{l - dl){xnix)die'^\x - y))Xm{y)u*{y)'F{y)\ (36) 

and estimate the various factors on the RHS. 
First, we use, for j < 4, 

\d^e'^\x-y)\<Ce-^"^ (37) 
on the support of Xm and Xn- To bound F{y), we need the bounds (for < i < 2) 

liaVlloo, ||aV||oo<Ce (38) 
which follow from (3), the Proposition and (20,27,28), and the bound 

\d'w{x)\ = \Xnix)d'wix)\ <Y,I dk\x^w{k)\ 

n n 

< ^(1 + n^)-^ I dk{l + e)-^\\w\\ < C\\w\\ (39) 
To extract the factor (1 + rn^)~^, we show that 

/ \x^mw{y)\dy < (40) 

and similarily for v* and u* (with i > for the latter) with \\w\\ replaced by e. 
Indeed, let e C^(R) be such that (t>rnXm — Xm- 

Then 



J \Xmd'w\ = J \Xm(pmd'w\ < (^J 



\xmdM') '<c[j Ix^^^r) ' < (41) 

For V* and u* (with z > for the latter) the bound holds with replaced by e due to 
the explicit expressions and bounds (3), (20) and the Proposition. 

Now we bound (36): the x integral is controlled by Xn{x) or its derivatives, the s 
integral is less than and we use (38) and (39) for u*{yy and all factors in F{y) except 
one, for which we use (40): 

Pnm < L2c'+^+^e-l"^-'^l(l + m")-^e'+^||«;||^^ (42) 
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with M the total number of factors of v* and its derivatives and of derivatives of u* in 
F, and N similarity for w. We have I + M + N >2. 

(B) Let now \m — n\ < 2. The difficulty is that we do not have (37) for s close to zero. 
But we do not have to control a sum over m G Z, and we can use Fourier transforms. 
Let us denote (pm.u*'' by where 0^ is as in (A). Then 

Xnd'am{k) = j dsdpdqxnik - p){ipye'''''' fm{p - q)XmF{q). (43) 

Let us consider the various factors on the RHS. Since x is with compact support, 
we have 

\Xn{k-p)\ = le-'^'-'Hik-p)] < Q{1 + \k-p\'r' (44) 

for any /. 

For frn: note that 

J 1(1 + {-d^Y)<PmU*'ix)\dx < CriCeY (45) 

for all r, whence 

\fm{p -q)\< a(C6)'(l + {p- qfT'- (46) 
Also, Jds\p\^e-'P < CL"^ if j < 2, so, provided we can show 

\x^F{q)\ < C^+^^e^||w||^(l + m^)-i(l + q')-' (47) 
we can perform the convolutions in (43) to get 

Pmn < L2c'+^+^(1 + m')-'e'^^'\\w\\^. (48) 

Using (42), (48), and the contraction mapping principle as in the proof of Theorem 1, we 
find a solution w{x,t) such that 



\w 



i;L')-e^^'-'^''gi.)\\<C,,aeie+\\g\\). (49) 



The first term comes from = 0, / + > 2 and the second from / + Af > 1, A > 1. To 
apply the contraction mapping principle, we used ||e*^*^^)'^ g\\L < C\\9\\ which is proven 
hke (52) below. 

To prove (47), use 

\x^'w\ < (1 + m^)"^(l + k'^y^\\w\\ 

and 

\d^{k)\^\J2x:d^w{k)\<^^ 

together with similar bounds for v* and the derivatives of u* and perform the convolu- 
tions as in equation (3.13). 
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(b) Contraction. Let us now turn to the RG (29), i.e. 

Rg = Lw{L-,L^). (50) 

Note that because of the derivative in K (see (30)), R preserves the condition ^(0) = 
i.e., by (29), ^ / w{x,t)dx — and so, J{Rg){x)dx — J g{x)dx — 0. We will show that 

ll^^ll<^ll^ll + CL,ae(6 +11^11). (51) 
Indeed, by (49), all we need is to show 

po^|| = ||L(e(^^-i)^^5)(L-)||<^|b||. (52) 
It is now easier to work in the x representation. We have 

Rog = j G{x,y)g{y)dy (53) 

with 

G{x, y) = (47r(l - L-2))-ie-2(i-^"')"'(^-^"'2')'. (54) 
We need to bound 

||i?o^|| = sup (l + n^)| / dxe'^'^ f dy{l - dl){xn{x)dl,G{x,y))g{y)\. (55) 

k,n,i<2 J J 

We consider first \n\ > ClogL. Then 

|...| / dxdy\{i-dl){xn{x)diG{x,y))xm{y)g{y)\ 

m '' 

< CEe-l"-^""^l(l + m^)-iy|| < CL-^TT^ (56) 
m \i. -\- n ) 

for L large enough. We used (41) (with d^w replaced by g) and (54). To control the 
sum, we used \n\ > ClogL for \n — L~^m\ > and (1 + m"^)"^ for \n — L~^m\ < 
which implies m > L-y-. 

For |n| < ClogL, we subtract G{x,0) from G{x,y) in (55). We may do this for free, 
since / g{y)dy — 0. Then, 

|...| <^ / dxdy\(l-dl)(xn{x)di(G(x,y)-G(x,0)))xm(y)9(y)\. (57) 

m 

For the terms with \m\ < L, we use 

\d^{G{x,y) -G{x,0))\ < CL-^e-l"l(l + |m|) (58) 
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and thus, using (41) again, 

y < y e-i'^i^L < cL-i^lL. (59) 

,^ ~ 1+ m3 - l + n4 ^ ^ 

\m\<L \m\<L ' i i 

For the terms with |m| > L, the subtraction in (57) is not necessary, (41) suffices. 

y <C y < CL-^\\g\\ < (60) 

\in\>L \'m\>L 

for L large, since \n\ < ClogL. Hence, (52) follows. 

(c) Iteration. To conclude the proof, we need to iterate R. Going back to (31), we study 
Rl,k„- From the expression (30) for K and from (32), we see that the only change to 
the n = 1 analysis above will be a change in (51): 

C 

\\RL,K„gn\\ < j\\9n\\ + ClML'""^ + \\9n\\)- (61) 

Indeed, all the terms in (30) are irrelevant, except some of those with / > 1, = 1, M = 
0, 1, which are bounded by CL,ae||5'ri||. Hence, we get (recall that \\g\\ < C\\f\\ by (28)) 

\\RL^,Kg\\<CL,aL-^'-'^^\\f\\ + e) (62) 

which allows us to conclude the proof as before. □ 
Proof of the Proposition. To solve (10), we set 

h = Atp (63) 

and solve 

h = dia*i(Pl + A~^h)di(Pl + ^'^h)) = ^0 + N{h) (64) 

in the Banach space (13). We need some information on the operator A. In fact, A is 
just the harmonic oscillator in disguise: 

e^Ae-'r = -d^ + — + -. (65) 
16 4 



Hence, the kernel of A ^ is readily computed from Mehler's formula |T^. The result is 



A'\x,y) = M{x,y)= i dtMt{x,y) (66) 

with 







_i -i, 

e 2 1 (x-e 2y)^ 

Mt{x,y) = ^ -re i— ' . (67) 

7r2 1 — e ' 2 
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In N{h) we encounter terms dA and d'^A . We write the latter as d'^A — 
{-A - lxd)A-^h ^-h- \xdA-^h. Thus, 

N{h) = R{A-^h, dA-^h, xdA-^h) (68) 

with R : ^ C analytic near zero. Thus all we now need is 

Lemma. The operators A"^ , dA"^ and xdA"^ are continous in the norm (13). 

Proof. We want to show 

1(1 + \x\)d\A-^h){x)\ < Ce-'r\\h\\j^ , i<N+l. (69) 

Let 

M = Mtdt + J Mtdt = M^^'^ + M^'^^ (70) 

For M^^) we have 

\d\M(^^h){x)\ <l^dtl dy\d'Mt{x,y)\e-'^\\h\\N (71) 
Prom (67) we have the estimate, for any 5 > 0, 

\d'Mt{x,y)\ < Cr^V5(l - e-*)-5(^+i)e-^*''T-e-'^ (72) 

Thus, using 

f l-S (x-e~^yf g-* I j 1-6 7? 

/g 4 i_e-t 8 = 20r((l - 5) Z7+-)~^e 4 l+(i_25)e-t (73^ 

J 1 — e * 2 

we see that for r large enough and b small, 

1(1 + \x\)d\M^''^h)[x)\ < C\\h\\Ne-^ (74) 
For M^^) we have (let i — N + 1, the worst case), integrating by parts, 

\d^+^{M^^)h){x)\ = I £ dte""^ J dyd^Mt{x,y)d''h\ (75) 

and thus by (72) and (73) this is bounded by 

|a^+i(M«/i)(x)| < [dtj^^^^^^e-'^^H^fit) (76) 

where f{t) = (^(1 - e"*) + (1 - 6)6-^^/'^ is 0(1) for all t. Now, using e"* < (1 + 
we have 

l-S x^ x^ x\ 1-26 , 



4 l + (l-25)e-*- 8 8'2-2(5 + t' 
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and so we get from the t integral 

— r+'v'^ — — 

yielding 

1(1 + \x\)d^+\M^^^h){x)\ < C||/i||7ve-^ . 

For j < N + 1, write 

|a^ (MW/i)(x)| = ij^dte^^ J dyMt{x,y)d^h\ (79) 
and get, from (72) for i = and from (78), 

1(1 + |x|)9^'(M«/i)(x)| < C||/i|Ue-'r. (80) 

The Lemma is proved. 

Solving (64) is now trivial. We have 

\\MN<Cae''+' (81) 

and 

||iv(/^)||^<c.(6'^||/.|U + 11/^11^) 

By the Lemma then 

\mN<c\\hu 

and the Proposition is proved. □ 

The perturbed Burgers' equation 

Let us now discuss the perturbed Burgers' equation, (3.36). To understand the rele- 
vant fixed point, consider first the case H = 0. Then, by the Cole-Hopf transformation 
0, the Burgers' equation reduces to the heat equation : set 

V^(x,t) = e^--"^"'*^'^ (82) 

Then, ijj = ip" if 

u = u" + {u^y (83) 

Now, lim^^_oo '^/'(a:, 1) = 1, lima;^+oo '(/'(x, 1) = e-^-°= '"'•^'^■''^^^ so that we look for a fixed 
point for ip of the form (3). Using the fact that the inverse transformation of (82) is 
u{x,t) = ^^^'l^ , we get the family of fixed points 



Ap'(r) d 
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for the transformation (2.4) with a — 1, where u solves (83). 

Note that the derivative in (84) explains why is a fixed point of (2.4) with a — 1 
while (/)q in (3) is a fixed point of (4) (where a — 0). 

Since e(— oo) = 0,e(+oo) = 1, 

/ + 00 
nix) dx (85) 



Now we can state 

Theorem 4 Let -ff : — > C &e analytic in a neighbourghood of with dp > and 
fix a 5 > 0. Then, there exists an e > such that, if \\f\\ < e, equation (3.36) with 
u{x, 1) = f{x) has a unique solution which satisfies 

hm t'-'\\ui-t'/\t)-t-'/'f*A-)\\=0 

t— »oo 

(the norms here are defined by (3.3) 

Proof. The proof is very similar to the one of Theorem 1, except for the change of fixed 
point. Let 



with 



so that, by (85), 



f{x) = fX{x) + go{x) (86) 



/ + 00 
fix) dx (87) 
-oo 



/+0O 
go{x) dx = ^o(O) = (88) 
-oo 



_j,2/4 

We have \\gQ\\ < C||/||: using e'{x) — ^ ^ i we see that f\^ and all its derivatives 
are integrable; using ||/i||oo < ||^||i we have 

ll/loll<^l^o| (89) 

Also, 

1^1 < Cll/ll (90) 

for 11/11 < e, because, by (87), |ylo| < C*!!/!!! < C||/||. For the last inequality, we 
use ll/lli < C||(l + |x|)/||2, by Schwartz' inequality, and ||/||2 + ||a;/||2 < C||/||, by 
Plancherel. 

The local existence is proven as before. To study i?/,, let us write the solution u{x^ t) 

as 

u{x, t) — u*^^{x, t) + v{x, t) 
24 



with 

Since u*y^^{x,t) solves (83) with u\^{x, 1) = fX^ix), v{x,t) satisfies : 

v^v"+ + vf - {uXfy + H{u, u") (91) 

with 

v{x,l) = go{x) (92) 

Now, 

Rl /(•) = Lu{L-, L') = /!„(•) + L viL; L'). 

Write it as 

RLf{-) = fW + 9i{-) (93) 

with gi{x) dx = 0. 

This means that we define Ai (see (85)) by 

1 _|_ /I, /•+00 

1 -r ^0 J—oo 

But, by (88), (92), Jv{x,l) dx = and, by (91), \fj v{x,t)dx\ < J \H{u,u' ,u")\dx < 
\\H{u^u',u")\\ which, by a bound similar to (3. 14), (3. 17), gives 

/+00 
v{x,L')dx\<CL,H\\f\\' (95) 
-oo 

So, inserting (95) in (94), we have, for ||/|| small, 

\A^-Ao\<Cr.,H\\^t (96) 

We also have 

ll^ill < CL''\\go\\ + C^,H\\ff < L-^'-'^\\f\\ (97) 

This is like (3.25) : writing (91) as an integral equation we have Lv{L-,LP') = i2o9'o+ 
rest; i?ofi'o contracts because of (3.4) and the rest is the sum of two terms which are 
bounded using (3.13): the first, coming from 2{u*j^^vy, is less than Ce||5fo|| < L^^H^foH, 
(because, using (90), Hm^qHl + ||(itlo)'||L < C\\f\\ < Ce), and the second (coming from 
{v^y + H) is bounded by CL,H\\ff- 

Now the iteration is exactly as before. The H term runs down with L""*^^ so that, 
using (95), we have 

\Ar,+^-Ar,\<CL,HL-'"'" 
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and, as in (97), 

\\gn\\<CL-^^-'y^\ 

The rest of the proof is as in Theorem 1. □ 

5. Discussion 



Now, we want to put our results in a more general framework, relate them to other 
works, and explain the analogy with the theory of critical phenomena. The discussion 
will be heuristic or based on previously known results, and we shall limit ourselves to 
the family of heat equations with absorption: 

u = u" - VJP (1) 

where p > 1 is not necessarily an integer. In our terminology, p > 3 is irrelevant, p = 3 
marginal, and p < 3 relevant; 3 has to be replaced by 1 + in dimensions. 

To understand the p < 3 case, observe that (1) is invariant under the scaling trans- 
formation (2.3) for a = Finding a fixed point for that transformation amounts to 

finding a solution of the form u{x,t) = t~~ f*(-^), where /* is a solution of the ODE 
(replacing ^ by x) : 



/" + Ixf + ^-fP = (2) 
2 p — 1 

The following results are known about the solutions of (2): 

1) For any 1 < p < 3, there exists an everywhere positive solution of (2), which 
has almost Gaussian decay at infinity P, [TT| . 

2 

2) For any p > 1, there exists a solution which decays at infinity like |x| p-i 
|TT|, |TB|]. Note that is integrable only for p < 3. 

Rather detailed results are known on the basin of attraction of the various fixed 
points; we state them loosely, see the references for more precise statements e.g. on the 
type of convergence; also, in each case, the decay in time of u{x,t) is t~ , as in (2.10), 
where the exponent a is related to the fixed point as in (2.3, 2.4). 

1) For p > 3, and initial data u{x, 1) non-negative and integrable, the asymptotic 
behaviour of the solution is governed by /g, as in Theorems 1, 2, with logarithmic 
corrections for p = 3 ||lT], |12|. This is a global result {u{x, 1) is not assumed to be small); 
our results are perturbative, and restricted to integer p (but a general non-linearity F), 
but do not require u{x, 1) to be pointwise positive and hold (when p > 3) also if one has 
in (1); in that case, smallness of u{x, 1) is necessary since large initial data blow up 
in a finite time iTDI, ISI, M • 



2) For p < 3 and u(x, 1) non-negative and having (suitable) Gaussian decay, the 



asymptotic behaviour is governed by the non-trivial fixed point /j" |Tl], [17 . 
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3) If one starts with u{x, 1) non-integrable and decaying at infinity like witli 
< a < 1, then, for p > 3, the relevant "Gaussian" fixed point is /* where f^ik) = 
|^|a-ig-A: ^ -^jiiQii^ for any a, is a fixed point of (2.4), where u solves the heat equation. 
This fixed point has the right decay at infinity. Now is relevant, marginal, or 

irrelevant according to whether p<l + ^, p=l + ^, oTp>l + ^. One knows |T^ 
that, for a non-negative initial data, the solution converges to /* for p > 1 + ^. For 
p = 1 + it converges to . For p < 1 + f , the solution converges to a solution constant 

in space, {p — l)^~t^~ , which solves (1) without the diffusive term: ii = —u^, and 
which can be viewed as a (somewhat degenerate) new fixed point. 

To see the analogy with the theory of critical phenomena, consider an Ising model 
or a 0^ theory, on an dimensional lattice, at the critical point. > 4 is like p > 3 
here: 0^ is irrelevant and the behaviour at the critical point is governed by the Gaussian 
fixed point. For N = 4, (p'^ becomes marginal, and the Gaussian behaviour is modi- 
fied, like here in Theorem 2, by logarithmic corrections. However, this is not true for 
every marginal perturbation. In 0^ theory, like here for p = 3, this happens because 
the marginal term becomes irrelevant when higher order terms are included: An and, 
therefore, /„ go to zero which is the same thing as having a coupling constant in front 
of the term going to zero. In particular, this higher order irrelevancy depends, like 
in (p^ theory, on the sign of the perturbation. For in (1), the solution blows up 
1^, In point 3) above, the marginal perturbation {p = 1 + -), leads to a non-trivial 
fixed point, /2 instead of /*. Also, in the Burgers' equation (3.36), the marginal term 
remains marginal to all orders and the solution is governed by a new fixed point, which 
is however easy to write down (see the end of Sect. 4). 

For A^ < 4, 0^ becomes relevant and one expects the critical behaviour to be governed 
by a non-trivial fixed point, whose existence is however much harder to establish than 
here. Another analogy with field theory concerns the constant A in (7): this is like a 
" renormalised" constant whose corresponding "bare" value is Aq = /(O). One of the 
problems encountered in proving (7), instead of just a bound on u{x, t), is that A, unlike 
Aq, may depend in a complicated way on / and F, and is not known explicitely. This 
usually limits the power of ordinary perturbation theory: one may try to expand u{x, t) 
around A^t'^f* which has the wrong constant and the perturbation series may not 
converge. One of the advantages of the RG method is that it allows to "build up" A 
through a convergent sequence of approximations, as we did in the proof. Also, note 
that in the marginal case, the renormalisation of A drives it to zero and produces the 
logarithmic correction in (39). 

Finally, one may also interpret in RG language the results on the existence of singular 
or very singular solutions of (1). A solution is singular if, as t — > 0, u{x,t) becomes 
concentrated on a point, and it is very singular if, moreover, / dxu{x,t) diverges when 
t — s> 0. For the family of equations (1), one knows that, for 1 < p < 3, there exists both 
a singular and a very singular solution 0, the singular solution is the fundamental 

one while u{x,t) = fi{^) is very singular, since fl is integrable. Besides, the 
fundamental solution behaves, for t — 0, like the corresponding solution of the heat 
equation (^~2/*(-^)). On the other hand, for p > 3, there does not exist any singular 
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solution 

This can be understood in RG language: start with some non-singular data at t = 1; 
a singular solution would be obtained by solving (1) backwards in time and letting t 0. 
This amounts to running backwards the RG flow; hence, the stability of the fixed points 
is inverted: /g becomes stable for 1 < p < 3 and unstable for p > 3. So, for p < 3, 
one would expect the solution to be attracted to /q when t —>■ 0, thus explaining the 
presence of a fundamental solution. On the other hand, for p > 3, one could only go, 
when t — *■ 0, towards /2 , which, because of its power law decay at infinity, cannot be a 
singular solution (it is not integrable for any t). This correspondence between t — »• and 
t — oo is similar to the correspondence between "ultraviolet" and "infrared" behaviour 
of fixed points in field theory and in critical phenomena. 
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